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Abstract 




Starting from a covariant cycle-averaged Lagrangian the relativistic oscil- 
lation center equation of motion of a point charge is deduced and analytical 
formulae for the ponderomotive force in a travelling wave of arbitrary strength 
are presented. It is further shown that the ponderomotive forces for trans- 
verse and longitudinal waves are different; in the latter, uphill acceleration 
can occur. In a standing wave there exists a threshold intensity above which, 
owing to transition to chaos, the secular motion can no longer be described 
by a regular ponderomotive force. 
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A so-called ponderomotive potential <3> p is induced by the oscillatory motion of a free 
charge. This potential plays a dominant role in atomic physics (e.g. multiphoton ionization) 
and laser plasma dynamics (e.g. parametric instabilities, self focusing, beat wave accelerator, 
fast ignitor U). It was shown independently by several authors || that for a monochromatic 
electromagnetic field of arbitrary space dependence, 

E(x,t) = $lE{x)e- lujt , (1) 

the oscillation center dynamics of a charge q is governed by the so-called ponderomotive 
force f p , 

Eq. (2) was obtained from a first order perturbation analysis of the Lorentz force [[J around 
the oscillation center and is therefore subject to the usual smallness constraints of certain 
parameters. In order to obtain some weak generalizations of this expression, a variety of 
different approaches was chosen [4-6]. Again, they are characterized by (i) a perturbation 
analysis of momentum, (ii) harmonic fields of type Eq.(l), and (iii) a small ratio of oscillation 
amplitude to wavelength A. Recently, however, in connection with the existence of new lasers 
capabable of delivering ultrahigh irradiance a real need for relativistic expressions for / , 
not bound by such limits, has arisen. One of the aims of this letter is to show under which 

conditions f n exists and what forms it assumes. 

p 

The relativistic Lagrangian L(x,v,t), v = dx/dt, of a charge q in an arbitrary electro- 
magnetic field E = — V$ — dA/dt is given by 

2 

L(x,v,t) = -— + qv- 7 = (1 - v 2 /c 2 )- 1/2 . (3) 

7 

When an oscillation center exists, the transformation to action-angle variables S = 
S(x, t), rj = rj(x, t) is possible (e.g. rj = k ■ x — ut in the case of a travelling monochromatic 
wave). The action S and the angle rj are both Lorentz- invariant. The motion of the particle 
is governed by Hamilton's principle, 

fV2 (If 

tS = 5j^ L(x(r ] ),v(r ] ),t(r ] ))-d V = 0. (4) 



From the Lorentz invariance of S and rj follows that the Lagrangian C(r]) = L(dri/dt)~ l 
is invariant with respect to a change of the inertial reference system. Assuming that 77 is 
normalized to 2ir for one full cycle or period of motion, the cycle-averaged Lagrangean £ ? 

Uv) = ^ £ +2w jC(v'W, (5) 

depending only on the secular (i.e. oscillation center) coordinates Xq, Vq through rj is defined. 
The oscillation center motion is governed by the Lagrange equations of motion, 

dt dVn dXn 



2 



with Lq = Codrj/dt. To demonstrate this assertion we prove the following 
Theorem: The validity of Eq.(4) implies 
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'ij 



(7) 



Thereby A" = (j]f — i]i)/27r is the number of cycles over which £ undergoes an essential 

change. The symbol o(A^ _1 ) means: "vanishes at least with order 1/iV". 

Proof: Let the variation be an arbitrary piecewise continuous function A(r/). The n- 

th cycle starts at r\ = rj n , at which for brevity we use the symbols A n = A(rj n ), 

dC /drj n = (dC /dr)) n=rjn . If the same quantities refer to an intermediate point r\ < r\ a < 

rj n + 2tt we write A a and dCo / di] a and omit the index n for the interval. In leading order 

holds: 
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/ C (rj + A)dr] - 27r£ (^n + A„) - / C (r])dr] - 27rC (rj 

JVr,. UVn 
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^A( V )d V - 2vr^A r 
vn orj n dr) n 



drja 
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In the last step the mean value theorem is used. The function A (77) is arbitrary. Therefore 
at i] = i] n A n = A a can be chosen now without affecting A a . With this substitution the 
leading order gives the result 



C dr] 



<(2vr) 2 E 
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drj. 
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In this last step it is essential that dCo/dr] is a smooth function (in contrast to A which 
is generally not). Now, N = min ^-\C 0max / (d£ /dr] n )\ is chosen, i.e. over N cycles £ c 
changes at most by Cq. It follows that 



max 



C dr] 



1 r 

< x max A 

~ N 11 



(8) 



Performing the variation of this inequality leads to Eq.(6) with the replaced by a function 
/ not larger than |jComax| x max |A|/iV 2 . • 
In order to understand what inequality (8) means let us specialize to a case of the 
averaged Lagrangian Co not depending explicitly on time. Then the Hamiltonian Hq = 
p ■ v — L (x , v ), where L = —Cq(uj — k ■ v ), owing to dH/dt = dH/dt = —dL /dt = 0, 
expresses the energy conservation H = E = const. A straightforward estimate shows that 
the uncertainty / in Eq.(6) leads to an energy uncertainty AH/H < 2tt/N. This means that 
Eq.(6) is adiabatically zero and the total cycle-averaged energy is an adiabatic invariant in 
the rigorous mathematical sense 0. For — ► 00 Eq.(6) becomes exact. Perturbative 
averaging of a Lagrangian was used in M . 



3 



The relativistic Hamiltonian of a point charge in the electromagnetic field follows from 
Eq.(3), 



H = p-v-L=[m 2 c A + c 2 (p-qA) 2 } + $, (9) 

with the canonical momentum p = dL/dv = 7777/u + qA. Its numerical value is the total 
energy E = / ymc 2 + $. Considering a monochromatic wave in vacuum we can set $ = 0. 
This motion is exactly solvable || . The cycle-averaged quiver energy in the oscillation center 
system is given by 

H$ — mc 2 

with a = 1 for circular and a = 2 for linear polarization. If the effective mass m e ff = 
—jC,o r Yo{dv/d't)/c 2 is introduced, Lq = Codr]/dt shows that in an arbitrary inertial frame in 
which the oscillation center moves at speed Vq, Lq and Hq are those of a free particle with 
space and time dependent mass m eff , 

Tfl 2 

L {x ,v ,t) = — , H (x ,p ,t) = 7o^effC 2 ; (11) 

7o 



.2 \ 1/2 



2\ -1/2 



^° = ( 1 ~ ^2 ] ' P° = 1om cS v . 



Expressions (11) hold in any electromagnetic field in vacuum in which an oscillation center 
can be defined. In the special case of Eq. (10) m eff = (1 + q 2 A ■ A / am 2 c 2 ) 1 / 2 . 

In the oscillation center system, i.e. in the inertial frame in which at the instant t v (t) = 
holds, the ponderomotive force follows from Eqs.(6) and (11), 



fN = dPo = dL 
p dr dxn 



-c 2 Vm eff . (12) 



The force f is given the index iV since it is a Newton force (co-moving system); r is the 
proper time. The Minkowski ponderomotive force F p , valid in any inertial system, is 

F P = (ft + ^{f N p ■ v )v , ^v ■ fA . (13) 
\ c / 

Hence, the three-dimensional relativistic ponderomotive force f at any oscillation center 
speed Vq is the Einstein force, 

f p = -- jv 7V m eff + -^(v ■ VV)t)o) • (14) 
To 1 v$ J 

In the non-relativistic limit Eq.(2) is easily recovered. 

To see the power of the Lagrangian formulation, Eq.(6), we calculate /_ in a non- 
relativistic Langmuir wave of the form E(x,t) = E(x, t) sin(fcx — cut) with slowly vary- 
ing amplitude E. In lowest order the potential is $(x, t) = (E/k)cos(kx — u)t) and 
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L = mv 2 /2 — g<3>. In the frame comoving with xq the particle sees the Doppler-shifted 
frequency Q = uj — kv ( plus higher harmonics which are not essential here). With the 
periodic excursion ((t) around x the potential is <3>(x, t) ~ $(x ,t) + ((t)d§/dx . From 
this L = ttiVq/2 — aE 2 /Q 2 , a = q 2 /4m, results in lowest order. With this Lagrangian it 
follows from Eq.(6) for E = E(x) (no explicit time dependence) that 

f dv (i-V )(i-3V ) d & 

/„ = m— — = —a —E . (15) 

P dt u 2 {l-V ) i -6aE 2 /mv 2 dx 

In the last expression Vq is the oscillation center velocity normalized to the phase velocity 
v v = u/k, Vq — Vo/Vtp. The formula shows that f p changes sign when the particle is injected 
into the Langmuir wave with a velocity Vq exceeding v v /3. In regions where the standard 
expression for f p , Eq.(2), always exhibits repulsion the more exact treatment can lead to 
attraction. Eq. (15) was also derived in a more formal but physically less transparent 
manner in ||. Uphill acceleration in an electron plasma wave of increasing amplitude E(x) 
is confirmed by the numerical solution of the exact equation of motion (Fig. la). If E(x) 
grows indefinitely the exact ponderomotive force changes sign again and the particle stops 
and is finally reflected; the corresponding path in phase space exhibits a hysteresis (Fig. lb). 



Uphill acceleration, occasionally observed in other context [fTOj] , is a general ponderomotive 
phenomenon. 

Strong gradients of amplitude and, as a consequence, large ponderomotive forces appear 
when two or more waves superpose. Of particular importance is the case of a standing, or 
partially standing wave. We have numerically studied the secular motion in the field given 
by the vector potential 

A(x, t) = Ae y (smr] + — sin r/_), r]± = kx±uit. (16) 

For field amplitudes £ < 0.1 (non-relativistic case), £ = \q\A/mc, splitting of the motion 
into a fast and a secular, i.e. ponderomotive one, makes sense. In Fig. 2a an electron orbit 
is followed over 110 laser cycles in the standing wave of strength £ = 0.1 : the motion is 
regular and its dependence on the phase of the wave field is weak (narrow band of orbits). 
The corresponding Nd laser intensity (co> Nd = 1.78 x 10 15 s -2 ) is J Nd = 1.21 x 10 16 Wcm" 2 . At 
£ = 0.25 (I Nd = 7.6 x 10 16 Wcm" 2 ) some electrons starting near the maximum A(x) = 2A 
are able to escape from the ponderomotive potential well; the secular motion through the 
well exhibits a strong phase dependence. At £ = 0.5 (/Nd = 3 x 10 17 Wcm -2 ) the motion 
becomes totally chaotic (see Fig. 2b): slight changes in the initial time produce totally 
different orbits, a clear signature of chaos. The chaotic motion has its origin in the Doppler 
effect since the moving electron "sees" one wave as blue- and the other one as red-shifted. 
No oscillation center exists under such circumstances and a secular time scale may only 
build up on the statistical average. 

Electrons in a dense plasma behave differently since they are coupled to the massive ions 
by an ambipolar electric field E s . In Fig. 3 the calculation of Fig. 2 is repeated for £ = 0.5 
in a plasma of which the density is assumed to be such as to keep the oscillation center 
in a steady state position. In Fig. 3a the single arc-like (no longer 8-shaped) orbits are 
shown. The inclusion of E s in the equation of motion is essential for producing the regular 
orbits and the correct ponderomotive forces (Fig. 3b). For comparison the bare dotted line 
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is the ponderomotive force from Eq.(2) without including the space charge field E s in f p . 
We estimate from numerical runs that beyond 7n<i — 10 18 Wcm -2 no regular ponderomotive 
force exists in the dense plasma either. 

We conclude as follows, (i) When an oscillation center of motion exists the invariant 
cycle-averaged Lagrangian describes the ponderomotive motion in arbitrarily strong fields, 
(ii) The ponderomotive force in a monochromatic travelling, locally plane wave of arbitrary 
strength in any reference system can be expressed analytically, (iii) In a longitudinal wave 
uphill acceleration and phase space-hysteresis may occur. Finally, (iv) superposition of 
different modes of the same frequency leads to a limiting intensity above which transition to 
chaos occurs. The influence of dissipation (radiation losses, friction) on the ponderomotive 
force is also well understood now, but will be treated in a separate paper. 

*) On leave from Rand Afrikaans University, Johannesburg, South Africa. 
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FIGURES 



FIG. 1. Ponderomotive motion of an electron (q = — e) in a travelling Langmuir wave. In (a) 
the velocity V = v/v,p, v ip = ui/k, is shown as a function of position X = kx in the longitudinal field 
£(1 — exp(— 0.1X)) sm(X — T), £ = eE/mujv v = 0.1, for the initial conditions Xq = 0, V$ = 0.6 
(trajectory a). The oscillation center motion from Eq. (2) (dotted line b) shows a slight decrease, 
whereas the secular acceleration by Eq. (15) (dotted line c) is in good agreement with that resulting 
from trajectory a. The lower curve d shows the harmonic electric field as "seen" by the moving 
electron. In (b) the normalized kinetic energy Vq/2 is shown as a function of $ p /m e ^ from Eq. 
(2) for the linearly increasing field 2.5 x 10 -3 £Xsin(X — T), £ = 1, and injection velocities 
Vq = 1/3 (lower curves) and Vq = 0.43 (upper curves). The motion shows a pronounced hysteresis 
in phase space owing to different acceleration in co- and counter-motion. 

FIG. 2. Regular (a) and stoachstic motion (b) over 110 laser cycles of an electron placed 
with oscillation center speed vq = at the position X = kx$ = 397r/40 in a plane monochromatic 
standing wave; ordinate Y = ky points into field direction. In (a) the maxima of the pondero- 
motive potential <5 p are located at X = 0, tt; at X = ir/2 & p = 0; normalized amplitudes of 
counterpropagating waves are £ = eA/mc = 0.1. The electron motion evolves in a narrow band 
between the potential maxima. Band broadning around X = tt/2 is due to Doppler effect. In 
(b) <!>p (dashed line) is 25 times stronger (£ = 0.5), the electron motion is erratic, sweeping over 
21 maxima of & p during 110 cycles. Note also the enormous transverse excursion of Ay = 0.9 
wavelengths (Ay = 5.8). 

FIG. 3. Electron orbits and ponderomotive potential in a standing plane, monochromatic 
wave in a uniform plasma, (a) Electron orbits around oscillation center tied to the ions by the 
induced ambipolar field; field distribution as in Fig. 2(a), field strength £ = 0.5 as in Fig. 2(b). 
In (b) the segments parallel to X = kx indicate the orbit widths in X-direction; the connecting 
dotted line indicates the relativistic ponderomotive force. For comparison, the bare dotted line is 
f p from Eq. (2); the field is given in dimensionless units of eE/muc. 
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